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1. INTRODUCTION 

The high reliability of induction machines (IM) is the main advantage which allows these machines 
to occupy an important place in industrial applications. A more reliable modeling of the physical phenomena 
governing the operation of IM, allows the development and implementation of linear [1] and 
non-linear [2]—[6] control strategies, allowing IM to have performances comparable or superior to those of 
DC machine [7]. Well known simplifying assumptions (isotropy of the machine, parameters invariance, 
neglected saturation [8]) lead to state models with constant parameters [9]. However, the stator and rotor 
inductances and resistances are variable due to saturation [10]—[13], and temperature rise [14]-[18]. Design 
of linear control laws based on these assumptions, are limited and are valid for a given fixed operating point. 
Nonlinear control strategies are developed based on the linear model without saturation [2], [17], [18], or on 
the model with saturation [3], [11], [19], [20], which allow to follow the speed and the flux references with 
accuracy depending on the complexity and the knowledge of the IM parameters. 

In this paper, starting from a nonlinear model of the rotor flux [16], we develop a nonlinear control 
law based on the backstepping technique (NLB), taking into acount the fact that the inductances (stator, rotor 
and magnetizing) of the IM machine depend on the magnetizing current. We use the rotor flux orientation 
(FOC)-which is the most used in literature-to deduce the control voltages. Robustness of this control law, 
with respect to heating, will be investigated and discussed in this article. For the symbols used, the reader can 
refer to Table 1. 
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Table 1. List of symbols 


Symbol Designation 
Usy, Usy Stator voltages in the rotor flux-oriented reference frame 
isx İsy Stator currents in the rotor flux-oriented reference frame 
lnk Rotor magnetizing current (A) in the rotor flux-oriented reference frame 
|[¥-| =Lmlimr| Rotor flux magnitude (Wb) 
(P, lref Rotor flux reference (Wb) 
L,(L,) stator (rotor) inductance (H) 
Lim Magnetizing inductance (H) 
R,(R,;) Stator (rotor) resistance (Q) 
Lro = Lm — Lr Rotor leakage inductance (H) 
ToT Rotor time-constant, modified rotor time-constant 
wy Electrical speed (rad/s) 
Wg General reference frame 
Wrref Speed reference (rad/s) 
Tem Th Electromagnetic torque (Nm), Load torque (Nm) 
2 
o=1- Lm Total leakage factor 
LsL; 
E DC bus voltage (V) 
NLB Nonlinear backstepping controller 
FO Flux observer 
IMNS Induction motor with saturated 
IMNS Induction motor not saturated 
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2. MATHEMATICAL MODEL OF IM TAKING INTO ACCOUNT SATURATION 

In the literature, several articles have studied the magnetic saturation in induction machines 
(IM) [3], [21]—[26]. The control laws based on the saturation of the IM machine are: linear controls with 
update of one or more parameters, or a nonlinear control with constant or variable parameters [27]. The 
dynamic model of the IM that takes into account saturation of the iron core in [24] will be used to develop a 
nonlinear control law for speed and flux. Starting from [24], the mathematical model of IM, which takes into 
account the saturation effect is described by the following state as in (1). 


di i2 
ae = —Cylsy + (w, + cT, (wg = wr) isy + Czimrx — (CT, = Az fi Te) Wg ~ 37 -0r ) imry — 2 i = 
mrx 
di te fitsx 
lsy ; p . * r . 
dt = Cyl sy 7 (w, + CT (Wg ~ wr)) İsx + C3 lmry + (CsT, n azı fıTr)@g Z csTpr ) lmrx 


yA 


ls 

y 

— C2 an: +fiusy 
S 


dimrx 
“dt T Azak sx ~ Az2lmrx + A327; (Wy ~ wr )imry 
dimr 
e BIRA AE ko * A 
dt Azalsy — A2almry — 0327 ;(Wg = wr )imrx 
dwr 3 y š . 
ai.” a330, + f3 (imrxisy om aoe — faTı (1) 
i i 
In the rotor frame reference, which rotates at speed Wg = Wmr = wr + a22 —~., the model (1) becomes : 
lmrx 
disx P è i2, * i2 
dt = —Cylsy + Orlsy + (a22 + c2) A + Cglinrx — C2 5 ate fitsx 
t lmrx mrx 2 
di, isi fia i 
y . n SX°SY . + 1%21 : sy 
dt = —(ay, = C2 )isy = Wrlsx T (a22 F c2) —— isa + C3linrq F @rlmrx — C27 F fisy 
d tt lmrx a22 lmrx 
lmrx n . 
: dt = 32 (isy aki imrx) 
W . . 
Ss Az3@, + fatmrxtsy ar fT, (2) 


So, the electromagnetic torque is given by the (3). 
Tem z IÈ İmrxisy (3) 


Coefficients that appear in (1)-(3) are defined by the (4). 


š Rs 1-0 , 1 + 1-0 , 1 by 1 3- 12; p š ” 
ay, = + ajz = az, =L a a = c. ajı + aj (AL — 
11 S gr Gre 2 T gup’ 021 spe 022 = Fe 033 mh PE 2P 7 fs rae 1 11 12( 
2 
2AL"), c, = Aj2AL", cz = G3, fı + aj2(AL — AL*), AL = L — Lm, AL’ = SAL (4) 
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T andL are called modified rotor time-constant and dynamic magnetizing inductance, respectively, and are 
defined as in (5) and (6). 


an (5) 
_ al _ : dLm 
= dimri em Y lirl aiar (6) 


These equations can be rewritten in a more condensed form as (7) and (8). 


disx 


ie E; + fitsx (7) 
disy 
ae Fy F fiusy (8) 
With: 
= : 3 iĝy ; id, _ : ‘ 
EF, = —Cylsy + Wrlsy a (a22 F c2); + Czlmrx 7 C25 Fy = — (a11 = Co)isy = Wrlsxy T 
mrx mrx 
P ;2 
lsxlsy . . faazı . Usy 
(Qz2 + C2) ine lsa + C3 lmry i a22 Wrlmrx — C2 Lie, (9) 


2.1. Modeling of the rotor flux of IM 

In general, saturation curve is a hysteresis form. In our model, only the first magnetization curve is 
considered. Based on the procedure described in [8], the magnetization curve of the IM is shown in Figure 1. 
A mathematical model is described by the (10) [24]. 


[P| = a(1 — eo Fler!) + ylinrl (10) 
From (10), an analytical expression ofL,, is obtained by (11). 


RA] a(1-e7Flimrl) 
Lm = Zt = 


(1) 


limrl limrl 


r 


IV |- Roror Flux (Wb) 
Lm, L - Inductances (H) 


14 16 18 20 


. 8 10 12 
ie! - magnetizing current (A) 


Figure 1. Typical rotor flux curve |¥r|, magnetizing inductance (Lm) and modified inductance (L) 


Notes: 

1) Some advantages over polynomial interpolation should be noted: 

— This model has three parameters to be identified/optimized 

— The three parameters are physically interpretable (have a physical meaning) 
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It can be used to analyse and predict the behaviour of the machine for magnetizing currents that exceed 
the nominal value 

Inductances can be deduced by a simple derivation. An interesting physical interpretation of the 
coefficients a, B and y can be made. Indeed: 


Pad Lm = a6 +y: This means that the tangent of the magnetization curve at the origin is equal to 
imr|>0 
a+y. This value represents the initial magnetic state of the iron core. 


ii lim Lm = y: Thus yrepresents an inductance. It can be interpreted as magnetizing inductance when 
imrl? +0 


the machine is completely saturated. 


2.2. Flux observer (FO) 


The flux is usually difficult or inaccessible to measure. Therefore, it is necessary to construct a flux 


observer (FO) that makes the rotor flux available for the synthesis of the proposed nonlinear backstepping 
(NLB) technique. Let’s consider the two state equations of the magnetizing current of (1), and replace 
imrxand imry currents by their estimatest,,;, et Îmry. Thus, we rewrite the two components of the 
magnetizing current as (12): 


di */, 4 a 
EA 422 (iy imrx +T, (Og Or Jinry ) 
‘ (12) 
dimy * A a 
ti . A A 
at a22 (isy imry ~T, (Og — @p imrx ) 


The resulting errors are (Îmry = İmrx — lmrx€timry = İmry — Îmry), and as in (13): 


di, es : 
ms a22 (in T, (Og -0r Jiny) 
5 (13) 
dimry * (7 ee 
ari =a (as +T, Og =O, Jimi) 
Consider the following Lyapunov candidate function as in (14). 
1 fy x 
Vobs = 2 (as + nry) (14) 
Its dynamics is given by (15). 
Vons = Ümrx- Imr T tmry: ling = -az (Pmr + Cay) (15) 


Since a3, > Othen the FO is globally asymptotically stable. 
In the rotor flux-oriented reference frame, the magnetizing current (|im;| = imrx) rotates at the speedwn,, 
which allows to write as in (16) and (17). 


dimrx 


dt =. 22 (isx = Îmrx) (16) 
i 
Wmr = Or + a22 mA (17) 


Then the rotor flux will be deduced from (10), (11) and (16), (17). Figure 2 shows the block diagram of the 
flux observer (FO): 


Figure 2. Block diagram of the flux observer (FO) 
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3. DESIGN AND ANALYSIS OF SPEED AND FLUX BACKSTEPPING CONTROLLER (NLB) 
We are interested in rotor speed and rotor flux magnitude control of the IM described by (2). Actual 
speed of the IM must follow its reference speed whatever the operating mode of the IM is. Synthesis of this 
nonlinear control law is based on the backstepping technique. 
1“step: We define the following tracking error as in (18). 
€y = Or — Wrref (18) 
Deriving e, with respect to time, we obtain as in (19). 


ey = Wy T Orref = —A330, + fatmrxtsy K faTı pri Orref (19) 


The term fzimrxisy is the electromagnetic torque of IM, so the choice of this term as a virtual control appears 
natural. Thus, we define a new error: 


2, = fz İmrxİsy =u (20) 


The dynamics of z4, is obtained by deriving (20) with respect to time. 


2 
._d *y - 3? yes Im df. . x 
Z1 =< (frimrxisy)- 01 i) Tm dt Imrxtsy + L, T (imrxisy) =o 


(21) 


* 2 1 t 7 ; ; sate : 
= f3| a22 ca | +1 (is = imrx) isy timrx-fy@O |+ S3 fiimrxtsy 7 Ì1 
Im Lr 


The magnetizing current of the rotorimry is not measurable because it is difficult to access, we replace it, in 
the control law, by its estimate ?,,,,,. and we rewrite Z, in a condensed form (22). 


Ż1 = H1 + AUsy (22) 


With: 


2 1 dimrx - a . A 
M1 =f (é - ~) AL + 1) me isy + Îmrx fy | — ùr a1 = fo fiîmrz 


Notes: 


d d 
The dynamics of the magnetizing and rotor inductances (Emand) are given by (23). 


dLm _ dLm dimrx _ L-Lm dimrx _ AL dimrx 
dt dimrx dt imrx dt imrx dt (23) 
dLy _ d(Lor+Lm) _ dlm 
dt dt ~~ dt 


Rotor inductance is no longer constant, it depends on the magnetizing inductance, thus on the 
magnetizing imrx current. Leakage Lor inductance is assumed to be constant, this assumption is common to that 
of [16] and [3]. But the rotor inductance depends on the magnetizing inductance: L, = Lor + Lm, which makes 
it depends on the magnetizing inductance. Consider the following Lyapunov candidate function as in (24). 


1 
V, => (e? + 22) (24) 
The dynamics of is V, given by (25). 


V, = e16, +242, = e(—a330, +u — fat, — Orres) + z(e +u + Ay Usy) (25) 


Int J Pow Elec & Dri Syst, Vol. 14, No. 1, March 2023: 60-76 


Int J Pow Elec & Dri Syst ISSN: 2088-8694 Oo 65 
To have global asymptotic stability, must V, verify the following condition by (26). 
VA = —k,e? = d,z? (26) 


Where k, and d, are any positive real design constants. By identifying (25) to (26), we obtain the (27). 
es = —A330, + V1 — fat, — Orref (27) 
—d;Zı = ey + Hı + a1Usy 


2"4step: the nonlinear controller must follow the flux reference. We denote Ê, the estimate of rotor flux ¥.. 
So, we define the following error as (28). 


ez =|" — Eler (28) 


Then we deriveez with respect to time, we obtain (29). 


Ê, dim 


L= =Y] 


è =2)¥, 7 


|¥,|-2|¥,| [¥,| =2 


ref 


ref ref $ r 


K (29) 


2 P 
P, -2| Pler ig 


PAES : 
= 2a79Lmimrxisx — 2422 ief 


We choose the term as 22717, îmryisx a second virtual command and define a new error Z, as (30). 
= 25 | 
22> 2a22Lmîmrxisx — U2 (30) 


The dynamics of z, is given by (31). 


. dlm; . 2. di, 2- dix. 
z2 =4a22Lm 4 imrxisy + 2a22Lmisx MEy 2a22Lmimrx See V2 
dt dt dt 
ok P z 2 ms Ar . 1 
=2a22L}m [a3 (2AL +Lin Jiz (isx —burx ) +P, fx) | + 2az2Lm |F r| usx — 02 3 ) 
We rewrite (31) in a more condensed form (32). 
Zp = H2 + QzUsx (32) 


With:w2 = 2az2Lm|a32(24L + Lm)isx(isx — İmrx) + AAE) — ùz; @ = 2az2Lm|®,|. 
We consider the new Lyapunov candidate function defined by (33). 


Vz == (e3 + 22) (33) 
Its derivative with respect to time is given by (34). 

V, = e262 + Z2Ż2 = €z (v = 2a22|P,|° = 21%-|rerl¥l er) + 22(€2 + Uz + A2Usx) (34) 
To have global asymptotic stability, V,must satisfy by (35) 

V, = —k,e3 — dzz3 (35) 


Where k, and d, are any positive real design constants. 
By identifying (34) to (35), we obtain the (36). 


ee = v — 2az2|%,|° — 21% re Bleep (36) 


—d2Z2 = ez + My + AzUsx 


According to (27) and (36) the control voltages u,,andu,,, satisfy as in (36) and (37). 
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—e1-d1Z1—-H41 

Usy = ay (37) 
—e2—d2Z2—-U2 

uoy = EEE (38) 


Now we need to calculate the derivatives of the virtual commands v,and v2. Using (27) and (36) we obtain: 


ù, = —kyé, + A330, + fat, + Orref 

= —033(433 — ky)w, + f3(433 — ky bmrxisy — fa( Cass — kı)t, + 1)ė, + kiOrref + Örref 
: : ala so 2 $ 
Vz = —k2ċ2 + 4a52|¥,| P,| + 21% er + 21% lrer”, or 


2 


= 2(2az2 — k2)|¥,||¥| + 21% [er 


+2 (ltor + hee) |P, lref 


4. SIMULATION RESULTS AND DISCUSSION 

The controller gains are selected as d4 = 1.10°,k, = 7.104, k = 7.104, d3 = 0.1.10?.The choice 
of the controller gain and the adaptive gain is not arbitrary. So, if this choice is correct, we get a high 
performance of the NLB controller. In a first test, rotor flux and speed references are set to 1 Wb and 100 
rad/s respectively. At 1 second, nominal load torque of 15 Nm is applied to IM. Figures 3(a)-(e) shows 
simulation results of this first test. 

It can be seen that rotor speed reaches its reference of 100 rad/s (Figure 3(a), in about /6ms with an 
absolute error of less than 0.08 rad/s in Figures 3(b)-(c)), before and after the application of the load torque. 
At the moment of application of 15 Nm load torque, the absolute error of speed does not exceed 1.8rad/s, 
before it returns to the range of 0.08 rad/s. Rotor flux -necessary for the synthesis of the control law- is 
estimated using the nonlinear flux observer (FO) depicted in Figure 2. Induction machine (IM) produces a 
rotor flux that follows well its reference (Figure 3(d). Before applying load torque, the absolute flux error is 
about 9 mWb, while after applying this torque, this error decreases to 3 mWb, but in all cases this error does 
not exceed 10 mWb (Figure 3(d)). 

Electromagnetic torque developed by the machine compensates the load torque Figure 4(a). In this 
test it is assumed that the friction is part of the load torque, and it is assumed to be known or estimated, and 
its first derivative exists. Figure 4(b) shows the error between the electromagnetic torque and the load torque. 
This error is zoomed in to see the effect of applying the load Figure 4(c), it shows that this error does not 
exceed 1.5 Nm before and after the load torque was applied. A torque peak is observed at time 1 second 
which is quickly cancelled by the control. i,, current component evolves in the same way as the rotor flux, 
and isy current component evolves in the same way as the electromagnetic torque. This result is common to 
the classical FOC, and it is also confirmed by in (3) and (10) as shown in Figures 3(d)-(e). 

Stator currents have a profile that confirms the torque and flux in (3) and (16). Indeed, in steady 
state, i „current component controls the rotor flux (in (2) and (10)) and Figure 4(d) shows this dependence. 
i,ycurrent component is responsible produces the electromagnetic torque (in (3)), because in steady 
state(imrx = is, )and therefore the electromagnetic torque depends essentially on isy as shown in Figure 4(e). 

As long as inductances L,,and Lare constant, currents i,,andi,,are also constant. Variation of rotor 
time-constant T, is also due to the variation of IM inductances (L,,L,,Lim,...). Figure 5(a) shows the 
variation of the rotor time-constant, the magnetizing and modified inductances obtained by the first test. 
Initially, both inductancesL,,, andL start from, Lm = L = 0.4695 H, which explains the peaks of current and 
electromagnetic torque observed at the beginning of the simulation, then they decrease until they stabilize, at 
no load (T, = 0), to the values :Lm = 0.1477 H and L = 0.028 Hrespectively. After applying load torque, 
the values of inductances increase slightly and become respectively:L,, = 0.1514 H and L = 0.030 H. 
Figure 5(b) gives the two time-constants introduced in the IMWS model. Initially the two time-constants 
described in the IMWS model, have the following values: T, = T,; = 0.3107 s, then they decrease to 
T, = 0.103 s and Tř = 0.0199 s. After applying load torque, these values slightly increase to 
T, = 0.105 s and Tř = 0.0212 s. From these results, it can be concluded that the model that takes into 
account the magnetic saturation, is valid to describe the variable behavior of inductances and the rotor time- 
constant (at constant rotor resistance). Figures 3 to 5 show, that rotor speed and flux are well controlled by 
the proposed nonlinear control law (NLB), despite the variation of several parameters such as the rotor (L,), 
stator (Ls) and magnetizing (Lm) inductances and the rotor time-constant (T,). 
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Figure 3. Simulation results for references | yj|,ep=1 Wb and ©, rep =100 rad/s; (a) rotor speed, (b) absolute 
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speed error (|w, — Wrref DE (c) zoom around the moment of application of the load torque (between 0.8 and 


1.2 seconds), (d) measured rotor flux of IM, and (e) absolute flux error (| l= [Se lee 


) 
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Figure 4. Simulation results for references | y;|,e=1 Wb and ©, ref =100 rad/s; (a) electromagnetic torque at no 
load and when applying a 15 N.m load torque, (b) the torque error (T, — Tem), (c) zoomed view of the torque 


error (T, — Tem), (d) and (e) stator currents i,,,and i,yrespectively 
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4.1. Influence of the rotor resistance variation 

In a second test, the speed of the machine is fixed at 100 rad/s and the rotor flux at 1 Wb. The 
objective of this second test is: 

— to compare the proposed control law in the case where the model of the induction machine takes into 
account saturation and in the case where the inductances are assumed constant; 

— to test the robustness of this control strategy (NLB) with respect to the heating of IM (increase of rotor 
resistance R,). Simulation results are shown in the Figure 6. 

Figure 6 shows the speed of the induction machine for different values of the rotor resistance, 
considering constant rotor speed and flux. It can be seen that the increase of the rotor resistance (by heating) 
has no influence on the steady state speed, neither in the IMWS case Figure 6(a) nor in the IMNS case Figure 
6(c). However, the influence of the rotor resistance is clearly seen in transient state. The increase of the rotor 
resistance allows the speed to reach its reference faster as shown in Figures 6(b)-(d). Indeed, if R, = 0.5R,, 
the motor speed reaches their reference at 0.02 seconds, and if R} = 2R,, , the speed reaches its reference 
after 0.01 seconds. This is justified by the fact that the rotor time constantT, decreases with increasing 
resistance, so the speed response becomes faster. 

Figure 7 shows the rotor flux produced (measured) by IM when taking into account the saturation 
(IMWS case) Figures 7(a) and 7(b) and when assuming no saturation (IMNS case) Figures 7(c)-(d) following 
a variation of the rotor resistance. In steady state, we note that by increasing the rotor resistance, the real flux 
produced by IM stabilizes at its reference of 1 Wb see Figures 7(a) and (c)) However, in transient mode, the 
rise time of each flux depends on the rotor resistance value. 

Indeed, forR, = 0.5R,n, the flux reaches the reference after 0.25 s in the IMWS case Figure 7 (b), 
and after Js in the IMNS case Figure 7(d). Whereas for R, = 2R,n, the flux reaches its reference more 
quickly, but it presents an overshoot. This overshoot increases with increasing rotor resistor value. In 
transient state, the difference between IMWS and IMNS cases is clear. In the IMWS case, the maximum 
overshoot is about 20 %, but in the IMNS case, this overshoot reaches 70 %. This is due to the fact that the 
magnetic saturation limits the overshoot to 1.2 Wb (flux of the saturated machine). However, in the IMNS 
case this limit does not exist because the flux is assumed to be linear (constant inductance). On the other 
hand, the flux response is faster in the case with saturation (IMWS) than in the case without saturation see 
Figures 7(b)-(d). 
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Figure 5. Motor’s parameters variation during the first test (a) magnetizing and modified inductance and 
(b) rotor and modified time-constants 
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Figure 6. Resulting speed for references |y,|re=1 Wb and ©, ref =100 rad/s; (a) IM speed in the IMWS case for 
increasing of rotor resistanceR,., (b) zoomed view on time interval [0-0.02 s] in IMWS case, (c) speed in the 
IMNS case, and (d) zoomed view on time interval [0-0.02 s] in the IMNS case 
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Figure 7. Resulting flux for references | y4|,.=1 Wb and ©, ref =100 rad/s, (a) actual rotor flux produced by IM 
in the IMWS case, (b) zoomed view on the transient state IMWS case), (c) actual rotor flux produced by IM 
in the IMNS case, and (d) zoomed view on the transient state (IMNS case) 
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Figure 8(a) shows the electromagnetic torque developed by the machine as a result of variations in 
the rotor resistance and the modified rotor time-constant used in the IMWS model. In steady state the 
electromagnetic torque follows the load torque. There is an overshoot at the beginning of operation due to the 
fact that initially the rotor time-constant is very large compared to its steady state value see Figure 8(b). 
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Figure 8. IM electromagnetic torque for references |y,|re=1 Wb and ©, ref =100 rad/s, (a) different values of 
rotor resistance, and (b) rotor and modified time constants in the IMWS case 


Table 2 gives a comparison between the IMWS model and IMNS model, using two performances 
indexes: the first (IAE: Integral of absolute error JAE = f fe |e(t)|dt), and the second (ITAE: Integral of time 
multiplied by absolute error ITAE = i. tle(t)|dt). Performance index (IAE, ITAE) is a parameter on the 
basis of which we can decide the system accuracy and sensitivity. Table 2 shows that the IAE errors of rotor 
speed, in the IMWS and IMNS cases, are very close and it decrease as the resistance R, increases. The ITAE 
error of rotor speed is lower in the IMWS case, compared to the IMNS case. Concerning the flux, we notice 
that the IAE error is better in the IMWS case than in the IMNS case, and it decreases as the resistance R, 
increases. This error remains constant after R, = 1.5R,,,. For the ITAE, its values are lower in the case with 
saturation, than in the case without saturation. Also, it remains constant after R, = 1.5R,». 

These results confirm that the model, which takes into account the magnetic saturation (i.e, IMS) 
allows to describe with more accuracy the phenomenon of magnetic saturation. The proposed nonlinear 
control law (NLB) allows minimizing the speed and flux errors in the presence of heating in IM. 
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Table 2. IAE and TIAE in the case of the IMWS and IMNS model 


IMWS Model IMNS model 
Wr — Orref kA = Fr lrep Wr — Orref P, z |F leer 

R,./Rrn IAE ITAE IAE ITAE IAE ITAE IAE ITAE 
0.5 1.174 0.035 0.064 0.011 1.133 0.087 0.149 0.041 
0.75 1.059 0.035 0.032 0.009 1.026 0.101 0.054 0.009 

1 0.991 0.034 0.021 0.008 0.961 0.095 0.013 0.001 

1.5 0.910 0.038 0.020 0.008 0.890 0.084 0.054 0.006 
1.75 0.886 0.041 0.020 0.008 0.863 0.084 0.066 0.006 

2 0.865 0.042 0.020 0.008 0.847 0.081 0.076 0.007 


4.2. Influence of the variation of the speed and rotor flux references 

In this section, we vary both the flux and rotor speed references, assuming constant rotor resistance 
equal to nominal value (R, = 1.0R,.,). Figure 9 shows simulation results for variable reference of speed 
(20 > 40 > 60 > 100rad/s), rotor flux (0.2 > 0.4 > 0.6 > 1Wb) and load torque (2 > 4 > 6 > 
15Nm). Rotor speed and rotor flux produced by IM are well controlled, even if their references change 
abruptly Figures 9(a) and 9(b). This shows that the proposed control law can control the speed and flux of IM 
for different operating points. On the contrary, a PID controller cannot control the machine in these different 
operating points without resizing it. The electromagnetic torque of the machine confirms the previous results 
(Figures 9(c)). It is noted that at each transition between two operating points, the torque developed has a 
large value to quickly stabilize to the new operating point. In steady state, since the magnetizing current imrx 
and the stator current iyare equal, this lets us write: 


Tem = fr isyl sy (39) 
|| = Lins isx (40) 


Therefore, i,, current component controls the rotor flux (see Figures 9(d)) and the electromagnetic 
torque depends on the two stator current components isy and isy. Appropriate control of {,,component makes 
the torque depend only on i,,component (see Figures 9(e)). However, in the synthesis of NLB, we did not 
take this particular condition into consideration, but on the contrary, the torque as described by (40) depends 
(in steady state) on the product of the two stator current components and on the parameter fy = f(imry,Lm)- 
Remark: 


We can synthesize a control law without using the rotor flux orientation condition: (w, = Wmr = 
wr, + a22 =) and still obtain the same results. So, to optimize the synthesis, we can use this condition of 
Mrx 
flux orientation from the beginning and avoid any unnecessary calculation. 


4.3. Simulation configuration 
In all simulations conducted in this paper, the MATLAB/Simulink environment has been configured 
as follows: 
— Solver: ode8 (Dormand-Prince); 
_ Fixed point with a sampling period: Ts= 10 us; 
— Switching frequency: fa= 5 kHz. 
Table 3 gives the nominal values of the asynchronous machine used in this paper: 


Table 3. Nominal values of IM 


Parameter Value 
Rated Power Pnom 2.2kW 
Rated Voltage Vnom 220 V 


Rated frequency fnom 50 Hz 
Rated Torque T em nom 14.9 Nm 


Pole pair sp 2 

inertia moment J m 0.0067 kg.m? 
E (DC BUS) 650 V 
Leet 0.012 H 

Ren 2.902 

By 1.550 
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Figure 9. Simulation results for variable speed and flux references, (a) rotor speed:20 > 40 > 60 > 
100rad_/s, (b) rotor flux:0.2 —> 0.4 > 0.6 > 1Wpb, (c) load torque, (d) i,,stator current component and (e) 
is stator current component 
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5. CONCLUSION 

In this paper, we have synthesized a nonlinear control law based on the "Backstepping"” technique 
(NLB). Assuming the nonlinear saturated/unsaturated IM models, and consequently variable magnetizing 
inductance, the NLB accurately performs control of the rotor speed and flux without the need of several 
control loops. The developed control law gives satisfactory performances in terms of speed, flux and torque, 
even if the machine parameters change: change of the rotor time-constant, variation of the rotor resistance 
(heating) and the magnetizing inductance (saturation). In the literature, there is not enough study of the 
behavior of the control law with respect to the variation of the rotor resistance. Therefore, in this work, we 
have seen that the proposed control law (NLB) allows to control the speed, the flux and the torque even with 
the variation of the rotor resistance. The robustness of the control, with respect to IM parameters variation, is 
extensively tested and compared to the model without saturation. 
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